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ABSTRACT 


Aims. We study the X-ray emission of low-mass main-sequence stars to derive a reliable general scaling law between 
coronal temperature and the level of X-ray activity. 

Methods. We collect ROSAT measurements of hardness ratios and X-ray luminosities for a large sample of stars to 
derive which stellar X-ray emission parameter is most closely correlated with coronal temperature. We calculate average 
coronal temperatures for a sample of 24 low-mass main-sequence stars with measured emission measure distributions 
(EMDs) collected from the literature. These EMDs are based on high-resolution X-ray spectra measured by XMM- 
Newton and Chandra. 

Results. We confirm that there is one universal scaling relation between coronal average temperature and surface X-ray 
flux, Ex, that applies to all low-mass main-sequence stars. We find that coronal temperature is related to Ex by 
©or = O.llEx' 26 , where T cor is in MK and Ex is in erg s -1 cm -2 . 


1. Introduction 


The X-ray emission properties of low-mass main-sequence 
stars are strongly dependent on stell ar mass and rota¬ 
tion (|Vilhul Il984l : iPizzolato et all 1 2 0031) . Although the ba¬ 
sic physical mechanisms responsible for heating coronae to 
X-ray emitting temperatures are not well understood, it is 
known empirically that coron al temperatures correlate well 


with levels of X-ray emission ( Vajanal 1983; Schriiver etjd 
1984: ISch mitt et ah 19901: Schmitt et al. 1995tlGudel et al 


1997; ISchmittl 11997 Tjellgscl©g©al. 2005b ITelleschi et al 
2007 1. For example, ITelleschi et al.l (|2005l) analysed high- 
resolution XMM-Newton spectra of six solar analogues with 
different levels of magnetic activity and showed that Lx ex 
^coi° 5± °' 25 j where T COI is the emission measure weighted av¬ 
erage coronal temperature (given by Eqn. [T|). 

An important question is which measure of X-ray emis¬ 
sion correlates best with coronal temperature for stars with 
a range of surface areas and spectral types; the three obvi¬ 
ous candidates are X-ray luminosity, Lx, X-ray luminosity 
normalised by the bolometric luminosity, li x = Lx /Lhni, 
and X-ray surface flux, Ex = Lx/(A-kR%). ISchmittl (1993) 
compared ROSAT hardness ratio measurements for a sam¬ 
ple of main-sequence F and G stars with similar measure¬ 
ments for a sample of main-sequence K and M stars and 
concluded that all low-mass main-sequence stars have the 
same relation between sur face flux and spectral hardness. 
Similarly. IPreibischl dl997t) analysed ROSAT spectra for a 
sample of stars and argued that Ex is a better measure than 
Lx, though they did not consider the parameter Rx- 

The answer to which of these parameters is best corre¬ 
lated with coronal temperature could contribute to our un¬ 
derstanding of the magnetic processes responsible for heat¬ 
ing coronal plasma to X-ray emitting temperatures and to 
our understanding of the coronal properties of saturated 
stars. For example, if Rx is the relevant parameter, then 
this would indicate that temperature scales somehow with 


how far a star is below the saturation threshold in rotation, 
given that saturation happens at a single mass-independent 
value of Rx- It would further imply that all saturated stars 
have similar coronal temperatures. On the other hand, if 
either Lx or Ex are the relevant parameters, then it would 
indicate that among stars that lie above the saturation 
threshold, lower mass stars have cooler coronae than their 
higher mass counterparts. The answer can also important 
for our understanding of stellar radiation in X-rays since 
stars with hotter corona will emit more photons at higher 
energies for a given total X-ray flux. This could be impor¬ 
tant for our understanding of the influences of stars on the 
upper atmospheres of planets, given that the influence that 
high-energy radiation has on a planet is highly wavelength 
dependent dTian et all [20051 ICecchi-Pestellini et al.1120091 : 
iChadnev et al.1 2015H . 


In this paper, we analyse coronal temperatures for a 
large sample of stars with different masses to show that Ex 
is indeed the best indicator of coronal temperature and to 
derive a general scaling law that can be used to predict T cor 
for all low-mass main-sequence stars. In Section 0 we anal¬ 
yse ROSAT measurements of hardness ratio for a sample 
of nearby X-ray emitting stars. In Section [3] we calculate 
Tcor for a sample of stars based on high-resolution XMM- 
Newton and Chandra spectra and show how it correlated 
with Ex- 


2. ROSAT Hardness Ratio Measurements 

To explore the differences between Lx, Ex, and Rx as in¬ 
dicators of coronal temperature, we collect a large sample 
of ROSAT X-ray o bservations from the NEXXUS database 
(ISchmitt fc Liefkell2004(P1 . NEXXUS is a database of X-ray 
measurements for nearby stars and provides a comprehen¬ 
sive compilation of ROSAT measurements, mostly from the 
ROSAT All-sky Survey (RASS). We collect X-ray luminosi- 
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Fig. 1. Plots showing how different measures of X-ray emission scale with hardness ratio for stars of differ ent masses. The h ardnes s 
ratios and levels of X-ray emission were measured by ROSAT and collected using the NEXXUS database (ISchmitt fe Liefkeil2004fl . 
The blue and red points represent stars with masses below and above 0.65 Mq respectively. The three quantities of interest are 
X-ray luminosity (left panel), X-ray surface flux (middle panel), and X-ray luminosity divided by the bolometric luminosity (right 
panel). In order to calculate Fx and Rx we estimate the surface areas and the bolometric luminosities by assuming R+ oc M 9 ' 8 
and Lboi oc M 8 ' 9 . 



Fig. 2. Best fit parameter, X, against the fit parameter c. For 
each value of c, we fit Eqn. [5] to the sample of stars given in 
Table. [l]and calculate X from Eqn. [3] The three vertical dashed 
lines show the value of c that correspond to the different measure 
of X-ray emission and the dotted vertical line shows our best fit 
value of c. 


ties and hardness ratios reported in the database (specif¬ 
ically the HR1 values) for all stars with masses between 
~0.2 Mq and ~1.1 Mq with Position Sensitive Proportional 
Counter (PSPC) measurements (mostly as part of RASS). 
Since hardness ratio is closely correlated with the temper¬ 
ature of the emitting plasma, they can be used to measur e 
coronal temperature (ISchmitt et, all IT 99, 51 : ISchmittl Il997h . 
The stellar masses are estimated by converting (B — V) 
colours from the NEXXUS database using a conver sion de¬ 
rived from the stellar evolution models of An et all (120071 1. 

The correlations between hardness ratio and the three 
measures of X-ray activity are shown in Fig. [TJ The red and 
blue points are for stars above and below 0.65 M 0 respec¬ 
tively (repeating this analysis with smaller mass bins leads 
to the same conclusion). Although there is overlap between 
the two mass bins in the Lx and i?x distributions, at a 
given temperature, higher mass stars are clearly likely to 


have higher values of Lx and lower values of i?x than lower 
mass stars. On the other hand, the distributions for the 
two mass bins in the Lx plot overlap excellently, providing 
a good visual indication that there is one mass-independent 
relation between Lx and coronal temperature. 


3. The Correlation Between X-ray Surface Flux 
and Coronal Temperature 

In this section, we show in a more comprehensive way that 
there exists one universal relation between coronal tempera¬ 
ture and Lx and derive a simple scaling law between the two 
quantities. For this purpose, we collect coronal temperature 
estimates from the literature based on multi-temperature 
emission measure distribution (EMD) fits to high-resolution 
Chandra and XMM-Newton spectra. Our sample contains 
24 stars in total with masses ranging from ~0.1 M 0 to 
~1.2 M 0 , including the Sun. For each EMD, we calculate 
the coronal average temperature by assuming that the log¬ 
arithm of the average temperature is the emission measure 
weighted average of the logarithm of the temperatures o f 
each component dGiidel et al.ll2007t iJohnstone et al.ll20 1 4l) , 
i.e. 


log Tcor 


Yn EM i ^g T i 

EM to t 


( 1 ) 


where the sums are over all measured components of the 
emission measure distribution with each component having 
an emission measure of EM, and a temperature of Tj, and 
EMtot = JAEM*. Stellar radii, Lx values, T COI values, and 
the references for our sample of stars are given in Table. |T] 
We make the assumption that the coronal temperature, 
stellar mass, and Lx can be related by 

LxM: = aT c b or , (2) 


where a, b, and c are parameters to be fit to the sample. 

The interesting parameter here is c since its best fit value 
is dependent on which measure of emission correlates best 
with coronal temperature. If Lx is the best parameter, then 
c ~ 0.0; if Lx is the best parameter, then c « —1.6; if L?x is 

1 The NEXXUS database can be found at 

http://www.hs.uni-hamburg.de/DE/For/Gal/Xgroup/nexxus/index.ht: 
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Fig. 3. Average coronal temperature against Fx (upper panel), Lx ( lower-left panel), and Rx ( lower-right panel), for a sample 
of low-mass main-sequence stars. As in Fig. [U blue and red represent stars with masses below and above 0.65 Mq respectively. 
The black line in the upper panel shows the relation T « 0.11 Fx' 26 , where T is in MK and Fx is in erg s _1 . The two data points 
connected by the dashed line show the Sun at cycle minimum and cycle maximum. The values and references for these stars are 
listed in Tabic |T| 


the best parameter, then c « —3.9. We define our goodness- 
of-fit parameter, X, for a given set of values for a, b, and c, 
as 

* = £ 7< [log {Lx,iM^) - log(a) - b log(T ,)] 2 , (3) 

i 

where the sum is over all stars given in Table. [T] and 7 , 
is a parameter that we use to weight the importance of 
certain stars in our fit. In order to make the set of low mass 
stars have the same importance in our fit as the larger set of 
higher mass stars, we assume 7 ,; = 1 for all stars with masses 
above 0.5 M e and 7 * = 4 for all stars with masses below 
0.5 M q . To find the value of c, we fit the free parameters in 
Eqn. [2] for each value of c between -5 and 1. In Fig. [21 we 
show the best fit value of X as a function of c. The best fit 


value of c is -1.45, which shows that coronal temperature 
correlates with Fx much better than with Lx and Rx- 
That Fx is better than Lx and Rx can be easily seen 
in Fig. [3] and is most clear from the two lowest mass stars 
in the sample, SCR 1845 and Proxima Centauri. Our best 
fit line in the upper panel of Fig. [3] is given by 

f cor « 0.11F°' 26 , (4) 

where T COI is in MK and Fx is in erg s ” 1 cm~ 2 . Our resu lt 
i s con sistent with the relation provided by iTelleschi et all 
(l2005| ) for just the solar analogues in the sample. Despite 
the large range of masses, all of our stars fit this relation 
excellently, providing confirmation of our conclusion that 
one universal scaling law exists between Fx and coronal 
temperature. 
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Name 

R* 

(R©) 

Fx 

(erg s -1 cm -2 ) 

Tcor 

(MK) 

Reference 

Solar Min 

1.00 

4.44 x 10 3 

0.97 

1 

a Cen A 

1.24 

2.46 x 10 4 

1.49 

2 

a Cen B 

0.84 

5.60 x 10 4 

1.76 

2 

Solar Max 

1.00 

7.73 x 10 4 

2.57 

1 

SCR 1845 

0.10 

2.30 x 10 5 

2.30 

3 

j3 Com 

1 .10 13 

2.47 x 10 5 

3.89 

4 

Prox Cen 

0.14 14 

3.36 x 10 5 

2.70 

5 

£ Boo B 

0.61 

4.13 x 10 5 

2.23 

6 

70 Oph A 

0.85 

4.24 x 10 5 

3.39 

7 

36 Oph B 

0.59 

4.31 x 10 5 

2.71 

7 

36 Oph A 

0.69 

4.35 x 10 5 

2.82 

7 

70 Oph B 

0.66 

4.65 x 10 5 

3.23 

7 

e Eri 

0.78 

5.65 x 10 5 

3.48 

7 

7 T 3 Ori 

1.32 15 

8.61 x 10 5 

4.31 

8 

X 1 Ori 

1 .02 13 

1.41 x 10 6 

4.37 

4 

k 1 Ceti 

0.93 13 

1.70 x 10 6 

4.57 

4 

£ Boo A 

0.83 

1.73 x 10 6 

4.37 

6 

tt 1 UMa 

0.96 13 

2.05 x 10 6 

4.47 

4 

YZ CMi 

0.36 

2.73 x 10 6 

5.79 

9 

AD Leo 

0.37 16 

7.04 x 10 6 

6.39 

10 

EV Lac 

0.36 17 

7.44 x 10 6 

6.78 

10 

AU Mic 

0.84 18 

1.28 x 10 7 

7.09 

11 

AB Dor 

1 .10 19 

1.39 x 10 7 

9.32 

12 

EK Dra 

0.91 13 

2.39 x 10 7 

9.12 

4 

47 Cas B 

1 .00 13 

4.04 x 10 7 

10.72 

4 


Table 1. Properties of our sample of stars shown in Fig. [jT] 
From left to right, the columns correspond to the stellar radii 
that we use in our calculations, the Fx values, the T cor values, 
and the references for the studies from which these values were 
derived. In most cases, we take the stellar radii from the studies 
that reported the X-ray luminosities and coronal temperatures; 
in cases where this was not possible, the superscripts on the 
R + values give the ref erences. The r e ferenc es are as follows: 
l.lPeres et al.l (l2000fl : 2,lR aassen et al.l (l2003fl: 3.lRobrade et al 


ll201Cfl: 4. Tellesclhet^l 


6 . 

8 . 

10 . 


Wood fc Linsky 


Wooc^&Laming 


^Robrade_& Schmitt] 

12. lLalitha et all (120131); 13 . 


_ —. _ Gride] 

(2009)~ 15. iBovaiian et ahl 1201 



i __GudeDa_ajl 
7. I Wood fe Linskvl 
9. Raassen et al 


11 . ___ 
n 


(2004) 
(2006) 

_ ( 2007 1 

iMagee et all (20031 


14. Demorv et al 


Favata et al. 


17. Huenemoerder_eta TT l2oic?> 7~ 18. iPlavchan et al 


19. IStrassmeiein2ob9l l. 


4. Conclusions 


It is very clear that the emission measure weighted aver¬ 
age coronal temperature scales very closely with X-ray sur¬ 
face flux. Our scaling law between Fx and coronal tem¬ 
perature can be used to accurately estimate the coro¬ 
nal average temperature of any low-mass main-sequence 
star when the value of Fx is known. Alternatively, Fx 
can be roughl y estimated bas e d on t he star’s ma s s and 
rotation rate dPizzolato et al.1 1 200. 'll : IWright et al.l l201lh 
iReiners et alJ [2014D . The scaling law that w e derive could 
be us eful as input into coronal models (e.g. iJardine et al.l 
120021 and stellar wind models that scale win d tempera¬ 
ture with coronal temperature ( Holzwarth_fc_J ar dinel [2007; 

I Johnstone et alJl2015aH Johnstone et alJ 2015bl) . 

For a star with a given mass and radius, the level 
of X-ray emission is determined primarily by its rota¬ 
tion rate, with quickly rotating stars emitting more X-rays 


than slowly rotating stars unt il a certain threshold where 
the X -ray emission saturates dPallavicini et alJl!981i [Vilhu 
I1984D . In the unsaturated regime, a star’s X-ray luminos¬ 
ity scales well with its r otati on period approximately as 
Lx oc n 2 (|Reiners et al.1 120141 ), which combined with the 
assumption that R* oc M 98 , implies that 

T cor oc M- 0 ' 42 ^' 52 . (5) 

A similar result can be obtained using the scaling laws of 
IWright et al.1 (|2011l) . Given that the saturation threshold 
appears to approximately be at a single mass-independent 
value of Rx , our result implies that in the saturated regime, 
low-mass stars have cooler coronae than high-mass stars. 
Assuming Tboi oc M 39 and R* oc M 98 implies that among 
saturated stars 


To, 


M 


0.6 


( 6 ) 


Assuming a sa turation Rossby number (= P TO t/T c ) of 0.13 
(IWright et al.lf20lD and a convective turnover time given 
by t c = 15.49 (AR/Mq)- 1 - 08 days (IReiners et al.ll2014h . the 
mass dependent satiation threshold is approximately given 

by 


Lip, 


= 13.53 



(7) 


where we use FIq = 2.67 x 10 -6 rad day -1 . 

The mass dependence of coronal temperature in the sat¬ 
urated regime is an interesting consequence of the fact that 
coronal temperature depends on Fx and not Rx and can 
have consequences for the influence of stellar high-energy 
radiation on the upper atmospheres of planets. Given that 
the location of the habitable zone is to first approxima¬ 
tion determined by Lboi, all planets within the habitable 
zones around saturated stars should be exposed to approx¬ 
imately the same stellar X-ray and EUV fluxes regardless 
of the central star’s mass. However, the spectra of high- 
mass saturated stars is likely to be harder (i.e. more pho¬ 
tons at higher energies) than the spectra of low-mass stars. 
The higher photon energies will mean that the radiation 
is likely to penetrate deeper into a planetary upper atmo¬ 
sphere, which could lead to differences in the atmospheric 
photochemistry and mass loss rates. 
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